Nonequilibrium mesoscopic conductance fluctuations 
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We investigate the amplitude of mesoscopic fluctuations of the differential conductance of a metal- 
lic wire at arbitrary bias voltage V . For non-interacting electrons, the variance {5g'^) increases with 
V. The asymptotic large-V behavior is {Sg"^) ~ V/Vc (where eVc = D jl? is the Thouless energy), 
in agreement with the earlier prediction by Larkin and Khmelnitskii. We find, however, that this 
asymptotics has a very small numerical prefactor and sets in at very large V jV^ only, which strongly 
complicates its experimental observation. This high-voltage behavior is preceded by a crossover 
regime, V jV^ < 30 , where the conductance variance increases by a factor ~ 3 as compared to its 
value in the regime of universal conductance fluctuations (i.e., at 1/ — > 0). We further analyze the ef- 
fect of dephasing due to the electron-electron scattering on (iJg^) at high voltages. With the Coulomb 
interaction taken into account, the amplitude of conductance fluctuations becomes a nonmonotonic 
function of V . Specifically, {hg ) drops as 1/V for voltages V 2> gVc , where g is the dimension- 
less conductance. In this regime, the conductance fluctuations are dominated by quantum-coherent 
regions of the wire adjacent to the reservoirs. 

PACS numbers: 73.23.-b, 73.63.Nm, 72.20.Ht 
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I. INTRODUCTION 

Universal conductance fluctuation g-'^i^i^ (UCF) repre- 
sent one of the most famous manifestations of quantum- 
coherent electron transport in mesoscopic conductors. 
Due to the quantum coherence of the electron motion 
over the entire sample, the variance {Sg'^) of the dimcn- 
sionless (measured in units of e^/ft.) conductance g is of 
order unity, with no dependence on the system size. The 
numerical coefficient depends on the spatial dimension- 
ality and on global symmetries of the Hamiltonian. In 
particular, for a quasi-one-dimensional system (a wire) 
and for preserved time-reversal and spin-rotation sym- 
metries one has ((5.g^) = 8/15 . 

In a pioneering paper** Larkin and Khmelnitskii (LK) 
predicted an enhancement of mesoscopic fluctuations of 
the differential conductance g = dl/dV at high bias 
voltages V. Specifically, they found that for V ^ Vc , 
where eVc = D/L/^ is the Thouless energy, D the diffu- 
sion constant, and the L the system size, the variance of 
the conductance increases linearly with V , {Sg"^) ~ V/Vc ■ 
This result implies, in particular, that regions of neg- 
ative differential resistance should appear in the IV- 
characteristics of a mesoscopic sample at sufficiently high 
voltages, V/Vc > g^- 

While early measurements^'^'^ did not confirm the pre- 
dicted enhancement of the conductance fluctuations by 
the bias voltage, more recent experimentsSiSiiSiii did ob- 
serve an increase of {5g'^) with V . This increase was in- 
terpreted as the linear behavior {Sg'^) ~ V/Vc predicted 
by LK. 

Inelastic scattering processes lead to dephasing, thus 
suppressing interference phenomena like mesoscopic con- 
ductance fluctuations. Indeed, it was found in all ex- 
periments that at sufficiently high voltages the ampli- 



tude of conductance fluctuations starts to decrease with 
V . However, no quantitative theory of the effect of 
Coulomb interaction on the conductance fluctuations has 
been developed. LK4 characterized the effect of inelastic 
scattering by a phenomenological inelastic out-scattering 
time Tin and the corresponding length Ljn = \DT\a ■ 
They obtained a suppression factor of {6g'^) of the form 
~ {Lin/Ly for the case Lin ^ L . No microscopic theory 
for Lin (which would produce, in particular, its explicit 
dependence on V) was presented in Ref. 0. The authors 
of Refs. [lOl and [ll| identified the ^-dependent dephas- 
ing length L^{V) as a relevant length scale induced by 
the inelastic scattering. They further argued that for the 
case of Coulomb interaction L^{V) can be obtained from 
the known resultiSii^ for L^{T) by a replacement T —^ V . 
They estimated the suppression factor of {Sg"^) by the in- 
elastic processes as {L^piV)/ LY, at variance with LK. It 
was found experimentally in Ref. Hol that {Sg"^) decreases 
as V~^ at large V, with 7 = 1.28 ± 0.12 . The authors ar- 
gued that the electron-electron scattering is not sufficient 
to convert the LK enhancement mechanism into such a 
fall-off, and ascribed this behavior to the electron-phonon 
scattering. 

The purpose of this paper is to develop a systematic 
theory of the voltage dependence of conductance fluc- 
tuations of mesoscopic wires. In Sec. ^ we calculate 
explicitly the dependence of {5g'^) on V/Vc for not too 
high voltages V <C gVc , when (as we will see in Sec. IIV|) 
the electron-electron scattering processes do not play 
any role. While we confirm the LK asymptotic result 
{5g'^) ~ V/Vc at y » Vc , we find that the corresponding 
contribution has a very small numerical prefactor. As a 
result, the LK asymptotic behavior sets in at very large 
V/Vc only and can hardly be detected experimentally. 
This high voltage behavior is preceded by a crossover 
regime, V/Vc ^ 30, where the conductance variance in- 



creases by a factor ~ 3 as compared to its value in 
the regime of universal conductance fluctuations (i.e., at 
V ^ 0). In contrast to the weak increase of the con- 
ductance fluctuations in the LK regime, this strong en- 
hancement of {Sg'^) in the crossover regime can easily be 
verified experimentally. In particular, we believe that it 
is this behavior that has been observed in the most recent 
and detailed experimental studyii2*ii 



the electron distribution function, 



In Sec, mil we consider the effect of finite temperature, 
T > eVc ■ We show that the temperature does not affect 
the conductance fluctuations at high voltages eV ^ T 
(up to a constant offset) but strongly suppresses them at 
lower voltages eV <^T . 

In Sec. IIVI we analyze the effect of dephasing due to 
the electron-electron interaction on {6g'^) at high volt- 
ages. We show that the Coulomb interaction makes the 
amplitude of conductance fluctuations a nonmonotonic 
function of V. Specifically, we find that {Sg"^) drops as 
1/V for voltages V ^ gVc , where g is the dimensionless 
conductance. Our results demonstrate that in the non- 
equilibrium situation the sample cannot be characterized 
by a single dephasing length, since the latter is position- 
dependent and strongly increases near the contacts with 
the reservoirs. Furthermore, we show that the conduc- 
tance fluctuations are dominated by quantum-coherent 
segments of the wire adjacent to the reservoirs. Accord- 
ing to our results, the maximal rms conductance fluctua- 
tions \/{Sg'^) are of the order of ^/g and occur at voltages 
of the order of gVc ■ Thus, the regions of negative dif- 
ferential resistance predicted by LK in the high-voltage 
regime do not appear. 

Sec. is devoted to a comparison with experimental 
data of Refs.lio andllll Our results are summarized in 
Sec. IVII Technical details of the calculations are pre- 
sented in several Appendixes. 



II. NON-INTERACTING ELECTRONS 
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Here D is the diffusion constant, v is the density of states, 
X is the coordinate along the wire, i?, A^ and K denote 
the retarded, advanced, and Keldysh Green's functions 
respectively, and the wiggly line denotes the external cur- 
rent vertex, with which an operator evx is associated. 
The derivation of Eq. (Q is outlined in Appendix^ The 
diagrams for current fluctuations can be obtained by con- 
necting two current diagrams by impurity ladders in all 
possible ways. We get the following six diagrams: 
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We consider a quasi-one-dimensional metallic sam- 
ple (a wire) of length L attached to two reservoirs 
with a voltage difference V. In order to calculate the 
conductance variance in the strongly non-equilibrium 
regime, i.e. at a large bias voltage V, we employ 
the Keldysh technique (see Ref. [Tj, for a review see 
Ref. llSf) . In the Keldysh formalism the ensemble 
averaged current can be represented with the help 
of the ensemble averaged diagonal Keldysh function 
{G^'ix)) = {Gf{x,x)) = 2mv[2n^{x) - 1] , where n, is 



c-diff: 
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As usual in the impurity diagram technique, the dia- 
grams consist of electronic vertices (Hikami boxes) con- 
nected by difFusons and Cooperons. The Hikami boxes 
take into account the possibility of inserting additional 
impurity lines connecting electronic Green functions of 
the same type (retarded with retarded, or advanced with 
advanced) without crossings. For the hexagonal vertices 
of the diagrams c-diff and c-coop, only the insertions con- 
tributing to the leading (zeroth) order in the external 
momenta are shown. Evaluating all the diagrams (the 
technical details of the calculation are presented in Ap- 
pendix^, we get the following result for the correlation 
function of mesoscopic fluctuations of the current^ 
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Here H is a rescaled diffusion propagator satisfying the 
equation 



dl 
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(l)i{x) - (f)2{x) }n^{x,x') = -S{x - x') 

(3) 

with the boundary conditions that H = if any of the co- 
ordinates x,x' is either or L, and 01^2 = ix/L)Vi^2 sue 
the potentials corresponding to the bias voltages Vi_2 • 



Equation jSJ is written for spinless electrons and un- 
broken time-reversal symmetry. Alternatively, it is valid 
for spinful electrons with a strong spin-orbit interaction, 
when only the spin-singlet channel contributes. If the 
spin rotation symmetry is preserved, Eq. (0) should be 
multiplied by a factor 4. If time-reversal symmetry is 
broken, Eq. (jSJ should be multiplied by an extra factor 
of 1/2 . It is worth emphasizing that, when the complete 
set of diagrams is taken into account, they combine in 
such a way that only the spatial derivative of the Keldysh 
function enters the result, Eq. Q. This is because only 
the electrons in the energy window V contribute to the 
transport. This property, which does not hold for indi- 
vidual diagrams but only for the sums of the diffuson or 
Cooperon diagrams, can serve as a useful check of the 
calculation. 

We assume first that the temperature of reservoirs is 
sufficiently low, T <^ eVc (or, equivalently, the thermal 
diffusion length Lx = {D/Ty^"^ is much larger than the 
system size), so that we can set T = . In other words, 
we will assume that at low voltage the system is in the 
regime of universal conductance fluctuations. The effect 
of a finite temperature, T > eVc , will be analyzed in 
Sec.lml 

When a quasi- Id wire of length L is connected to 
two perfect reservoirs with different Fermi energies, the 
electron distribution function in the wire will not be a 
Fermi function, but have a double-step shape. In the 
Keldysh formalism, this follows immediately from the ki- 
netic equation for the average Keldysh function, 
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with the boundary conditions 



(Gf (.)) = -2... X { ^ _V/(Yi% , ::' , (5) 

where / is the Fermi function. This equation implies that 
the Keldysh function is a linear function of the coordi- 
nate, 

(G^ix)) - -2^i^{l - 2/(e) + 2f [/(e) - /(e - eV)] } . 

(6) 
This yields (including a factor 4 to account for the spin) 

Vl/Ve V2/Vc 

{5IiVi)SI{V2)) = 16K' / dzi [ dz2 ^.,-z, , (7) 
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and 
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Here we have introduced dimensionless variables 
y = x/L , z = uj/eVc , and a ^ {Vi - V2)/K • 

To get the fluctuations of the differential conductance, 
we differentiate Eq. {T)) with respect to both voltages and 
set Vi = V2 = V. The result has the form (we set T = 
as discussed above) 



{6g6g) = {SgSg)^^ + {SgSg)^ + {SgSg}.^ 
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Let us now discuss the physical meaning of the three 
contributions (|11|) , (|12|l , and \i'A\ to the variance of the 
differential conductance. Consider first the low-voltage 
(linear) regime V ^Vc- In this case, the only effect of 
a voltage increment dV is to change the chemical poten- 
tial in the reservoirs, and thus to make the states with 
energies in the window dV available for electron trans- 
mission. The corresponding conductance fluctuations are 
given by the term Hll|l , reproducing the well-known UCF 
result,^ {Sg"^) = 8/15 . The two other terms vanish at low 
voltages and thus describe truly non-equilibrium effects. 
Indeed, out of equilibrium the effect of the applied volt- 
age is not only to change the chemical potential in the 
reservoirs, but also to alter the distribution of electric 
field in the sample. The shift dV of the voltage induces 
the variation d(j){x) = {x/L)dV of the electrostatic po- 
tential, and this variation affects all the electrons in the 
energy window V contributing to the current. This is 
the origin of the third term, Eq. H13|l . The second one, 
Eq. p2|) . is the cross-term due to correlations between 
the above two random contributions to dl /dV. 

The following comment is in order here. When consid- 
ering the conductance in terms of the scattering-theory 
formalism, one often replaces the problem by that of 
non-interacting electrons with different chemical poten- 
tials (or, equivalently, different concentrations) in the 
reservoirs. While the contribution {dgSg)Q is correctly 
reproduced in this way, the other two terms would be 
lost. What is not taken into account by such a non- 
interacting picture is that, for a system size much larger 
than the screening length, the main effect of changing the 
(electro-)chemical potential is not in altering the concen- 
tration but rather in changing the electrostatic potential 
in the reservoirs and thus modifying the potential profile 
in the wire. 

To evaluate the derivatives in Eqs. H12I) and H13|l . we 
expand 11 up to second order in a, 

(14) 



where HI = {~dy — i^) , and use the representation 
diagonalizing Ili , 
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where \n) = \/2sin(n7ry) and A„ = {mr)'^ — iz . The ma- 
trix elements of y in this representation are given by 
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We first find the asymptotic behavior of {5g'^) in the 
limit V ':^ Vc ■ It is easy to see that the leading contri- 
bution is of the form V/Vc and is governed by Eq. H13|l . 
Furthermore, only the term proportional to |np should 
be retained in S2, Eq. ©, to find the asymptotics. 
Changing the variables zi, Z2 to Z = [zi -\- Z2) / \2, ^ 
z = (zi — Z2)/V2 , using the expansion (|14|l . and eval- 
uating the integrals via the residue theorem, we get 
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where 



ci 



= 64 > - 

771, n>0 

= 7.785-10" 
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Therefore, although the large-F asymptotics is indeed of 
the form V/Vc obtained by LK, the corresponding nu- 
merical coefficient ci (not evaluated in Ref. |j) is very 
small. 

We now turn to the evaluation of the full crossover be- 
havior of {Sg"^) as a function of V/Vc ■ It is more cumber- 
some, but can be done in a similar way; see Appendix lO 
for details. The result is shown in Fig. ^ The striking 
feature of this plot is that the large-y asymptotics (|17f) 
becomes applicable only at very large values of V/Vc ■ 
The emergence of the linear contribution with a coeffi- 
cient given by Eq. ((T7|l can be seen for V/Vc ^ 50 but it 
becomes dominant only at V/Vc ^ 1000 . For V/Vc ^ 50 
the conductance fluctuations are well described by the 
asymptotic formula which includes, in addition to the 
linear term (|17|l . a ]/- independent contribution (which is 
formally of the next order but numerically dominates up 
to V/Vc - 1000), 
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with Co = 0.8964. On the other hand, this very slow 
increase of {Sg'^) (barely distinguishable from a satura- 
tion) at large voltages is preceded by a broad crossover 
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FIG. 1: Fluctuations (^g^) of the differential conductance 
as a function of the bias voltage normalized to the Thou- 
less energy, V/Vc . The inset shows the emergence of the 
linear behavior at very large V jVa ■ The dashed line rep- 
resents the asymptotic behavior, {Sg'') = 8/15 + co + ci V/Vc 
with Co = 0.8964 and Ci = 7.785 ■ 10"''. 



regime at V jVc S 50 , where the conductance increases 
by roughly a factor of 3. In the experiments, the increase 
of {Sg'^) is observed in the range V/Vc < 200 . As is ob- 
vious from Fig. ^ this enhancement should be mainly 
attributed to the crossover regime rather than to the LK 
linear asymptotics. We will return to the comparison 
with experiment later, after having considered the effects 
of finite temperature and of the interaction-induced de- 
phasing. 



III. FINITE TEMPERATURE 

In this Section, we address the effect of finite temper- 
ature on the non-equilibrium conductance fluctuations. 
We assume that the temperature is not too high, so that 
at low voltages the inelastic scattering can be neglected 
(the dephasing length L^{T) is large compared to the 
system size). The temperature however influences the re- 
sult via the change of the electron distribution function 
in the reservoirs. At a finite temperature, Eqs. ((TT|l ~ lfT^ 
are modified in the following way (here we switch back 
to dimensionful voltages for clarity). 
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f{ei + eVi) (/ is the Fermi func- 

Again the derivatives are taken 

at Vi^V2 = V ."" 

Evaluation of the full crossover (i.e. the conduc- 
tance fluctuations at arbitrary V/Vc and T/eVc) is too 
cumbersome. We thus restrict ourselves to the limit 
T ^ eVc , when the temperature strongly affects the con- 
ductance fluctuations at low voltages (in the opposite 
limit, T <^ eVc , the temperature is irrelevant for all volt- 
ages, so that the results of Sec. HTl apply*). This condition 
is reasonably fulfilled in most of relevant experiments, 
where the temperature is usually several times higher 
than the Thouless energy. 

We consider now the two limits of low {V ^ Vc) 
and high {V ^ Vc) bias voltages. In the first case, 
V/Vc —^ , only the voltage- independent contribution 
{dgSg)^ survives. For T ^ eVc the thermal smearing 
strongly suppresses this contribution compared to its 
zero-temperature (UCF) value. 
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In the opposite limit, when eV,T ^ eVc (the ratio 
eV/T can be arbitrary), we find that the high- voltage 
behavior H19|) is modified by the temperature in the fol- 
lowing way: 
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where the functions a (^) and ^ (^) have the form 
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It is worth emphasizing that for any T at sufficiently 
high bias voltages the T — result is recovered up to the 
missing offset of 8/15 . 

Adding the term {SgSg)f), Eq. ^^, to ^^, we get 
an interpolation formula which is parametrically justi- 
fied in the regime V ^ Vc and has the correct limit at 
V/Vc — > 0. This formula can thus be used as a conve- 
nient approximation for {Sg 5g) in the full range of volt- 
ages at T » eVc ■ The resulting voltage dependence of 
the conductance fluctuations is shown in Fig. |2] for the 
temperatures T = 20 eVc and T = 50 eVc . We compare 
these theoretical results with experimental data in Sec.lVl 
below. 
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FIG. 2: The voltage dependence of the conductance fluc- 
tuations as given by the sum of Eqs. 1231 and I24II . for the 
temperatures T = 20 eVc and T = 50 eVc . The zero-T re- 
sult of Fig. Q is also shown for comparison. It is clearly seen 
that the coefBcient of the asymptotic linear behavior is not 
changed by the temperature. 



IV. DEPHASING BY ELECTRON-ELECTRON 
SCATTERING AT HIGH BIAS VOLTAGES 

Up to now, we have neglected the dephasing induced 
by the inelastic scattering processes, which is justified 
for not too high bias voltages V. With increasing V, 
the electron-electron scattering becomes stronger, and 
eventually leads to suppression of the conductance fluc- 
tuations. The corresponding theory is presented in this 
Section. As in Sec. IIIII we will assume that the temper- 
ature at the reservoirs is sufficiently low (more precisely, 
T <^ geVc), so that the corresponding dephasing length 
is much larger than the system size, L^{T) ^ L . In this 
situation, the dephasing will be entirely due to a high 
applied voltage. 

In quasi-one-dimensional diffusive systems, dephasing 
is dominated by scattering processes with a small energy 
transferiSii^ oj <^T , which allows one to replace the dy- 
namically screened Coulomb interaction by a fluctuating 
classical potential. The same arguments apply in the 
presence of a high bias voltage, when the fluctuations of 
the electric field are characterized by the effective tem- 
perature Tes defined below. In view of the semiclassical 
character of the problem, we can use the Boltzmann- 
Langevin approach (see Refs. [l3 and [l8|) to determine 
the correlation function of the fluctuating field. The fluc- 
tuation of the local current density j can be expressed in 
the following way. 



^3 = -iqD5p + (TSE+ Sf 



(27) 



where the first term is the diffusion flow produced by 
the density fluctuation dp, the second one is due to 



the fluctuation SE of electric field {a = e^vD is the 
conductivity), and the last term represents extraneous 
(Langevin) current fiuctuations induced by elastic scat- 
tering. At equilibrium, these correlations arc found from 
the fiuctuation-dissipation theorem.^"' Out of equilib- 
rium, they have been studied in the shot noise context^ 
We only need the fiuctuations of a;-components of the 
current (those along the wire), 

{53'^J'\x,t)53'^^\x' ,t')) ^2cjT,^5{x-x')5{t-t'), (28) 

with the effective temperature T^six) given by 



Toff = / den^{x) [1 



si^)] 



(29) 



where ^^(x) is the (in general, nonequilibrium) distri- 
bution function. This function itself is strongly voltage 
dependent. 

For not too high V, the distribution function has a 
double-step shape lO; for higher voltages the inelastic 
scattering processes lead to thermalization of electrons, 
and the the distribution function assumes the Fermi 
shape with a local temperatureJ^ The boundary between 
the two regimes is determined by the energy relaxation 
time Tg: no thermalization takes place for eVc ^ r^^. It 
is important that the dephasing is dominated by small en- 
ergy transfers e <C Tcs , which implies^'' that the dephas- 
ing time T0 is much shorter than the energy relaxation 
time Te. Thus, even when the system enters, with increas- 
ing voltage, the regime of strong dephasing, eVc ^ t7^ , 
the energy relaxation still remains weak, eVc 3> t~^, in 
a parametrically broad range of voltages. From now on, 
we consider this regime, and find for the effective tem- 
perature 



Teff(x) = eV 



!('-!) 



(30) 



This effective temperature is position-dependent, taking 
its maximum in the middle of the wire and approach- 
ing zero at the contacts with the reservoirs. Using the 
continuity equation, 



Sp= - Sj , 



(31) 



and the relation between fluctuations of the electric fleld 
E = —V(f) and the density. 



(,) = M^ [/(«)(,) 



(32) 



where L/'°^ (jr-r'l) — e^/|j — r'\ is the bare Coulomb in- 
teraction, one finds 
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-iuj + Dq^ [l + j/C/(o)(g)] 



Sfx 



(33) 



Since for Coulomb interaction in metals t/C/'-*'-' (g) ^ 1 , 
Eq. I2ZI) reduces to 



aSE^{q) = -Sf^>^\q), q^O. 



(34) 



The q = -component of SE^ = —dx4> is zero due to the 
boundary conditions: there are no potential fluctuations 
at the perfect reservoirs. Therefore, the electric field fluc- 
tuations in the coordinate space are given by 



where j/< = min(2/, y') and y> — niax(j/, y') 



SE^ix) ^-- 



'5jr(^)-;^ I dx5jT\x) 



(35) 



Using Eq. (|3()|l . we obtain the correlator of the electric 
fields, 



{6Ex{x) 5Ex{x')) = - 



T^s{x) 5{x ~x') - jTcs{x) 



1 , 1 /■ 

— Tcff(x') -t- —^ dxTcsix) 



.(36) 



The correlator of the potential can be found by spa- 
tial integration. We introduce dimcnsionless variables 
y = x/L , T = eVct , 9 = 4>/eVc , and find 



2V 

(%,r)%',T')) = :T^<5(r-T')y<(l-y>) 

3Vc.9 



y{l-y)+y'{l-y') 



{y> - y<) 



,(37) 



As can be seen from (|37ll . the dimcnsionless param- 
eter V/Vcg is a measure for the importance of inelastic 
processes. Therefore, in the presence of the Coulomb 
interaction, the results of Section ITU hold under the con- 
dition V < Vcg ■ This condition excludes, in particular, 
voltages of order V ~ Vcg^ which are required to reach 
the regime where the fluctuating differential resistance 
would become negative in the absence of dcphasing. 



As we noted in the previous section, the part ~ jllp 
of Eq. I|13|) dominates at high voltages, so th at we can 
neglect Re 11^ . Following the approach of R,efs. ll3l2(i we 
express 11 as a path integral. 



iii'ri)=yi 



?2(T2) = yi 



Tr\n,\^ ^ J dy I dy2 J dT,dT2e''^^'-^-^ J V^Ah) J V^h) x 







«i(0)=y2 



?2(0)=»/2 



Tl 



exp ■^ / dti 





dU 



4 Vc 



Ci(ii)-z0(ei(ti),ti) 



6 , ^{Vi-V2, 



Vc 



6(<2)+«e(6(t2),t2) 



(38) 



We now perform the averaging over the random fields by close to each other, \£,i{t) — ^2(i)| ^ 1 • We thus expand 
using Eq. (|37|) . When dephasing is strong, conductance the terms in the action representing the dephasing in 
fluctuations are dominated by pairs of paths which stay \^i{t) ~ ^2(01 , which yields 



Trln, 



5l(Ti)=yi 



?2(T2)=yi 



= Jdyidy2 JdndT2e''^^'-^'^ J V^i{h) J V^t^) 



«i(0)=y2 




Ci «(^i - V2) 
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Ci(ii) 



«2(0)=iy2 



dt2 



4 Vc 



6(i2) 



2V 
Vc9 



Fiy2)\Ti - T2I 



2V 
Vc9 



inin(Ti ,T2) 

dt\Ut)-^2(t)\y2il~y2) 




(39) 



8 



F{y)^ly'{l~yf 



(40) 



Note that the factors ^(2/2) and 2/2(1 — 2/2) in the de- 
phasing terms could be equally well written as F{yi) and 
2/i(l— j/i) respectively, since we assumed that the distance 
between j/i and y2 is much smaller than their distance to 
the reservoirs, \yi — 2/2I <C j/i, (1 — j/i) . We will have to 
check the consistency of this assumption afterwards. 

The last two terms in the exponent of Eq. (|39|l rep- 
resent the dephasing effects. The first of them is the 
dephasing during a small time difference |ri — T2I be- 
tween the time extensions of the both diffuson trajecto- 
ries, while the second one is the dephasing during the 
motion of both diffusons. Since after the z-integration 
the characteristic values of |ti — T2I are of order Vc/V, 
the first term gives a contribution of order l/g to the 
exponent and thus can be neglected. 

Employing the Fourier transform. 



dre'^^S(T) , / dzE, = 27rE{0), (41) 



we can write the important part of Eq. H13|) (the part 
containing |np) for V ^ gVc as 



{^a^a)-. 



■yi-KVVc 



92 



dVxdV^ 



S(0) 



Vi=V^ = V 



(42) 



where S(0) is twice Eq. (|39|l with t\ set to ti . 

By changing to new variables Ci = (fi + ^2)/\/2 , 
(^2 = (fi — i,2)l \r^ we can perform the path integral over 
Ci , which is not affected by the random fields and just 
yields e(T)/\/2. Then 



{bg6g)^ = i2V2T:VVc 



d^ 



Wi^Va 



dy dTlyiO,T) 



Vl,2=V 

(43) 



where Iy(C,,T) obeys the equation 

= 5{05{t). (44) 



This can be made dimensionless by choosing new 
variables t = p^^^r , rj = p^^^x , X = p^^^^I , and 
V = i^/2{Vi - V2)/Vcp , where p=^/2y{l~y) 2V/Vcg . 
Performing the time integration over i, we get 



{h5g)2 = 87r 



4/3 



''"'i?«('" 



dy 



^='>iW-y)] 



7/3' 

(45) 
where Q satisfies (— 9j^ + |7/| + vif) Q{i]) — (5(r/) . We see 
that the resulting integral diverges at the upper and lower 
limits corresponding to the vicinity of the reservoirs. The 
reason for this divergence is as follows: Near the reser- 
voirs the expansion of the action in \(,i{t) — £,2{t)\ , used 



to derive Eq. H39|) . ceases to be valid, since the charac- 
teristic values of |Ci(t) — ^2(01 determined by the local 
value of the dephasing length L^(ii) become of order of y 
or (1 — y) . The above divergence shows that these small 
vicinities of reservoirs (where the effective temperature 
approaches zero and the dephasing length is larger than 
in the rest of the wire) in fact dominate the conductance 
fluctuations. 

The exact calculation thus requires using the full form 
of the correlation function H37|l . which makes the problem 
much more cumbersome, leading to a two-dimensional 
Schrodinger-type differential equation rather than to the 
one-dimensional equation 144|) . We can get, however, the 
result, up to a numerical prefactor of order unity, in a 
simpler way. Indeed, as follows from the above consid- 
eration, the integral must be cut off at a distance Ly^. 
to the reservoirs which is of the order of the local phase- 
breaking length, yielding the self-consistency condition 



Lye ^ L^{yc) = y DT^[Tcs{yc)] 



(46) 



where Tcs{y) = eVy{l - 
This yields a cutoff 



yc 



y) and t4T,s] ^ (Di^yr^g) 



1/3 



Vc9 
V 



1/4 



(47) 



Therefore, the amplitude of conductance fluctuations, 
which are dominated by the coherent segments within 
the distance Lye to the contacts, is given by 



{^9^9)2 



2Vc 

' V 



V 



(48) 



Thus, the variance of the conductance shows a maximum 
at V/Vc ^ g and then decays according to the 1/V power 
law. 

Strictly speaking, in order to show rigorously that the 
asymptotic suppression at high bias voltages is given 
by this power law, we still have to check that {SgSg)-^ 
and {Sgdg)^ remain smaller than {Sgdg)2 with increas- 
ing voltage. This turns out to be the case, as is shown in 
Appendix IdI 

A natural question to be asked at this point is whether 
the result (|48|l can be obtained from simple qualitative 
arguments. Indeed, we know that for the case of the 
linear-response transport, when the dephasing length L^ 
is governed by the temperature L^ = L^{T) , the con- 
ductance fluctuations of a wire in the high-temperature 
regime, L^(T) <^ L , can be estimated in the following 
way. The system is split into segments of the length L^ . 
The electrons within each of the segments remain phase 
coherent, yielding fluctuations of the conductance of the 
segment, {5g'^{L^)) - D/LlT, see Eq. (gg. A fluctua- 
tion Sg{L^) of the conductance of a segment induces a 
fluctuation in the total conductance 



Sg r^ (L^/L) Sg{L^) 



(49) 



Adding incoherently fluctuations of all L/L^ segments, 
one gets 



{5g')^{L,/L)'{Sg\L,)) 



L^T L 



(50) 



This agrees with the diagrammatic calculation^ (using 
Lfj, as a phenomenological parameter), as well as with 
the path integral calculationf22i yielding the dephasing 
length L^{T) induced by the electron-electron scattering, 

L^T) ~ (i?2i./r)i/3. 

A generalization of this consideration on the strongly 
non-equilibrium setup was attempted in Ref. |lO| but 
the obtained result (mentioned in Sec. was different 
from Eq. H48|) . The arguments presented in Ref. IQ, 
miss the following two non-trivial properties of the non- 
equilibrium dephasing, which distinguish it from its 
equilibrium, finite-T counterpart. First, as we have 
shown above, the dephasing length in non-equilibrium is 
strongly position-dependent, and the conductance fluc- 
tuations are dominated by the segments of the wire ad- 
jacent to the reservoirs. Therefore, representing the wire 
as an incoherent series of L/L^{eV) equal segments of a 
length L^{eV) is not justified. Second, the conductance 
fluctuations in a segment of a length L^ cannot be ob- 
tained simply using the results for a coherent sample of 
a length L^ subject to a voltage V{L^) — VL^/L . This 
is because the non-equilibrium distribution function in 
the segment L^ is not determined by a single parameter 
V{L^) but rather is governed by two parameters, V and 

We are now going to show that when these peculiarities 
are properly taken into account, the qualitative consider- 
ation does reproduce the result H48|) . Consider a coherent 
segment of a length L^ (to be specified later). Gener- 
alizing the consideration in Sec. IIIII the variance of its 
differential conductance can be estimated as 



{Sg'iL^)) 



1 



eVc{L4 



■Jde[AU,{e) 



(51) 



where eVc(i0) = D/L\ is the Thouless energy corre- 
sponding to the segment i^, and A/L^(e) is the differ- 
ence between the distribution function at the left and 
right boundaries of this segment. It is easy to see that 
A/i,^(e) is equal to L^/L in the energy window of the 
width eV , and is zero outside this window. This yields 



('55'(i0)) 



1 



L„ 



V 



e^.(.,)-nfj -kIti ^ <^^' 



According to Eq. H49fl . this induces fluctuations of the 
wire conductance of the magnitude 



(^5^> - ( ^ ) f 



(53) 



Finally, using that L^ — Lye with y^ given by Eq. (|47ll . 
we get {Sg^) ^ g^Vc/V ^ thus reproducing the result H48|) 
of the path-integral calculation. 



V. COMPARISON WITH EXPERIMENT 

In this Section we compare our results with the mea- 
surements presented in Refs. Hol and fill There, magne- 
toconductanc traces of a gold wire with a length L = 
1.5 /im and a Thouless energy of eVc = D/LP' = 3.4 /xeV 
were taken at T = 300 mK over a voltage range up to 

V ~ 3.7 mV. This yields T/eVc — 7.6 and the maximum 
value of V/Vc ~ 1000 . The dimensionless conductance of 
the sample was g k, 1400 . 

Fig. O shows our results for both the non-interacting 
limit (not too high voltages) and the strong-dephasing 
limit (very high voltages), compared with the measured 
data. It is seen that there is a reasonably good agreement 
between the theory and the experiment in the overall 
shape of {5g 5g){V/Vc) and in the magnitude of its en- 
hancement at the maximum. 

The main deviation is that the initial increase of the 
conductance variance is considerably less steep than on 
the theoretical curve. This discrepancy can be partly at- 
tributed to the following factors. First, we used for the 
full curve in Fig. O the formula H24|) which takes into ac- 
count the crossover at eV ~ T but not the crossover at 

V ^ Vc- Though the temperature was larger than eVc 
by a factor ~ 8 in the experiment, we know from Sec. ^ 
that the V/Vc crossover region is rather broad, extending 
up to V/Vc "^ 50 . Thus the V/Vc crossover is expected 
to overlap with the eV/T crossover under the experimen- 
tal conditions, making the total crossover region broader. 
Second, we did not take into account the effect of dephas- 
ing in the low-voltage region, where we used the result for 
noninteracting electrons. It is easy to see from Eq. H39|) 
that this will lead to a relative correction of order V/Vcg , 
which will provide a smooth matching with the strong- 
dephasing (high voltage) regime (^HJ- 

We also note that the linear behavior (|17|l could not be 
accessed reliably in this experiment, since it requires, in 
view of the very small value of the numerical coefficient 
(|18|l . very large voltages V/Vc ^ 1000 , while the dephas- 
ing was setting in at voltages several times smaller. A 
reliable observation of this linear behavior predicted by 
LK thus requires samples with still larger values of the 
dimensionless conductance g. 

At high voltages, the experimentally observed decay of 
the conductance fluctuations is consistent with our pre- 
diction of the 1/V decay, Eq. (|1S|) . see Fig.|31 Indeed, the 
value 7 — 1.28 ± 0.12 of the exponent of the power-law 
decay, {Sg^) ex V~^ is in good agreement with our result, 
7=1. We thus disagree with the conclusion of the au- 
thors of Ref. [IQi who argued, based on a naive estimate 
of the dephasing effect, that the electron-electron scat- 
tering is not sufficient to convert the LK enhancement 
mechanism into the observed fall-off, and ascribed this 
behavior to electron-phonon scattering. Experiments on 
samples with lower conductances might give more data 
in the strong-dephasing range and thus allow a more ac- 
curate experimental determination of the exponent of the 
power law. 
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FIG. 3: Compar ison o f our results with experimental data 
taken from Refs. llJlll The diamonds are the data from 
Fig. 5.4 of Ref.mi The solid line represents the interpolation 
formula obtained by adding Eqs. I23II and 12411 . multiplied 
by 1/8 to account for spin-orbit interaction and broken time- 
reversal symmetry. The dashed line represents the asymptotic 
suppression, Eq. 1481 . with the numerical coefficient chosen 
from the fit to the data. 



VI. SUMMARY 



In this article, we investigated the amplitude of meso- 
scopic fluctuations of the differential conductance of a 
metallic wire at arbitrary bias voltage V. For non- 
interacting electrons, we have found, by using the 
Keldysh diagram technique, that the variance {Sg'^) in- 
creases monotonously with V. The asymptotic large- 
ly behavior is {Sg"^) ~ V/Vc (where eVc — D/L^ is the 
Thouless energy), in agreement with the prediction by 
Larkin and Khmelnitskii. We find, however, that this 
asymptotics has a very small numerical prefactor and sets 
in at very large V/Vc only, so that it is very difficult to 
observe it reliably in the experiment. This high- voltage 
behavior is preceded by a crossover regime, V/Vc ^ 30 , 
where the conductance variance increases by a factor ~ 3 
as compared to its value in the regime of universal con- 
ductance fluctuations (i.e, at V^O). We further an- 
alyze, in the framework of the path-integral technique, 
the effect of dephasing due to the electron-electron in- 
teraction on {Sg'^) at high voltages. With the Coulomb 
interaction taken into account, the amplitude of conduc- 
tance fluctuations becomes a nonmonotonic function of 
V. Speciflcally, {Sg'^) shows a maximum at V/Vc ~ 9 
and drops as 1/V for higher voltages. In this regime, 
the conductance fluctuations are dominated by quantum- 
coherent regions of the wire of a length ~ L (gVc/VY^^ 
adjacent to the reservoirs. Our results are in good agree- 
ment with available experimental data. 
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APPENDIX A: KELDYSH DIAGRAM FOR THE 
CURRENT 

Here we give the derivation of Eq. JQ, which relates 
the current to the impurity averaged diagonal Keldysh 
function. 

We begin by expressing the averaged Keldysh func- 
tion (^G^ {xi,X2)) using the averaged diagonal Keldysh 
function (G^{x,x)). The off-diagonal components of the 
kinetic equations for the matrix Green's function add up 
to 



{e-H)G^^+G'^{€-H) = - [(2n, - \)G^ + G^(2n, - 1)] 

(Al) 
(here impurity averaging is implied). Using the gradient 
expansion, ^'^ one finds 

G^(x,p) = '-[ [2n,(x) - 1] G«(p)G^(p) 

+ id,n,{x) [G«(p)apG^(p) - G-^(p)apG^(p)]}. 

(A2) 

This is equivalent to the relation 

(G^(xi,a;2))- 
^ ^dx{G^{xi,x)){G'^{x,x)){G'^{x,x2)) .(A3) 



2'KVT 



This relation can be easily understood diagramatically, 
see Fig. 0] Using the relation^^ between the function 



K 



{G''ix,x)) 



K 



{G^{xux-2)) 



>5 



Xl 



X2 



FIG. 4: Diagrammatic representation of the relation between 
the averaged Keldysh function (G'^(a;i, 3:2)) and the averaged 
diagonal Keldysh function (G^(a;,x)), Eq. lASII . which is 
used in the derivation of Eq. (Q . 
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G< = {G'^ + G^ - G^)/2 and the electron density ma- 
trix p, 



G^{xi, t; X2,t) = iNp{xi ,X2,t) 
{N is the particle number), one finds 
(-e) f de d 



(A4) 



L{x) 



/£s<^'(^'-''L™- (*^' 



The contribution from the term G'^ — G^ is proportional 
to the spectral density and does not depend on the state 
of the system. Therefore only G^ contributes to the 
current |i^ 



(lA^)) 



i-e) d 



de 



{G^{x,x')) 



2m dx J 2ti 
Inserting Eq. (|A3(I into Eq. I|A6|I gives Eq. (QJ. 



(A6) 



APPENDIX B: EVALUATION OF THE KELDYSH 
DIAGRAMS 

In this Appendix we calculate the diagrams which arise 
in the Keldysh technique and contribute to the correla- 
tion function of currents (see SecE}. We will show that 
the sum of all the diagrams yields Eq. Q . 

We begin by presenting the expressions for the vertex 
factors (Hikami boxes). Because of their local character 
(the electron Green's function decays exponentially on 
the scale of the mean free path I, which is much smaller 
than the system size L), they can be calculated in the 
momentum space. 

The right Hikami box (containing two Keldysh ver- 
tices) of the diagram a-diff is of second order in the mo- 
menta: 

54 



2t:vDt'^ (-2^294 + ql + qf) 




where the qi are the incoming momenta and 
<Zi + 92 + 93 + 94 = ■ The left box (containing two 
velocity vertices) is of zeroth order: 



At:vDt^ 



Here D is the diffusion constant, v is the density of 
states and t~^ is the scattering rate characterizing the 
disorder strength. For the diagram b-diff the vertex fac- 
tors are of first order in the momenta: 

.92 .^ 93 





U 



ArnvDr^qi . 



91 ^ 94 

The vertex factors of a-coop and b-coop are the same 
as the vertex factors of b-diff up to extra signs arising 
from the direction of the Greens functions. 

The vertex factors of the diagrams c-diff and c-coop 
are only needed to zeroth order. Since these diagrams 
contain only one diffuson/Cooperon, this gives a result 
of the same order in the momenta as the other diagrams 
(containing one more diffuson/Cooperon) evaluated up 
to second order in the momenta. 



<3 



= 0. 



'Attj^Dt'^ 



Using these vertex factors and changing to real space 
representation, the expressions corresponding to the dia- 
grams are 



('5/«)a-diff 



{5ISI)^_ 



{sisi),_ 



diff 






diff 



dxidx2 

L2 



q2 a2 



dx\ 



V,,^,, V,,^,, G«(xi) G«(x2) -V,,-,, TTT^.,-., G«(xi) G«(x2 



dxi 



f— y . 

\47ri/ V 2'KVT 



^) {^^i^Dr^f j de,de2 j 



dxidx2 
L2 



d 



d 



a^^^^-^^ 9^^^^-^^ ^" ^''^^ ^^^^""'^ + '•'• 



, 



(Bl) 



(B2) 
(B3) 
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<^^^^>— P = {£)' \2..r 



('^^«)b-coop 



f— v> 

xAniJ \ 2'KVT 



' ^' (47r.i?r^)^ /ai&. /^^ (^^V.,..^ T',.., G^(xO G^C^^), (B5) 



<^^^^>^-™- - (4^)" V2... 



1-) i^in^^Dr^) I de^de, I 



dxidx2 



K I 



L2 



<5(xi - X2) T'ei-e, Gf^(xi) Gf^(a;2) , (B6) 



where 7'ej_e2(xi,X2) is the diffusion propagator satisfying the equation 

ei + (j)i{xi) - €2 - 4>2{xi) > Vei~e2{xi,X2) = - {ini^T^) S{xi - X2) , 



and 01.2 are the electrostatic potentials corresponding to the voltages V\,2 ■ Using the identity 



1 92 



4 9a:i9a:;2 



r'^ixi,X2) 



d^ 



9a::i9a;2 



7'(a;i,X2) - 



1 



AmyDr^ 



d{xi - a;2) 



'P{X1,X2) , 



we obtain 



(««>a-diff + b-diff = - {Dr'f j de,de2 



(««>a-.oon = -(^^')' /&ld62 



{5ISI)^_ 



I a— coop 



-coop + c— coop 



dx\dx2 



.,, , ^-Re7',,_,, 



-^G^{x{) -^G^(x2) 
oxi ax2 



L2 ■^-^^l^a^^"(-)a^^"(-) 



{DT^f J deide2 J 



dxidx2 1 -r^ ^2 



i9a;i 

I 



(B7) 

(B8) 

(B9) 
(BIO) 
(Bll) 

(B12) 



Adding up these equations and rescaling the propagator 
by a factor 2'kvDt'^ , 11 — 2tivDt'^V , we obtain Eq. |5)|, 

which is written in units of (e2/(27r)) = (e^//i) • 



there are contributions to the variance of the conductance 
which do not grow asymptotically as V/Vc but dominate 
in the intermediate regime. First of all, there is a con- 
stant contribution 



APPENDIX C: CROSSOVER FROM THE 

LINEAR RESPONSE TO THE 

LARKIN-KHMELNITSKH REGIME 



{S9Sg)o 



15 ' 



(CI) 



In this Appendix we calculate the variance {Sg'^) as a 
function of the bias voltage in the full range from the 
UCF regime {V/Vc <C 1) to the LK asymptotic regime 
{V/Vc 3> 1). In addition to the asymptotics Eq. (fT7|l . 



which gives the familiar UCF result in the limit of zero 
bias voltage. Second, there is a contribution from the 
term {dg Sg) -^ containing one energy integration, which 
can be evaluated as 



{^9Sg)i 



V/Va 

32 fdz^—Ti 
da 



n 



Ren^ 



a=0 



720 + 16 (^) - 360 (-1)^/' y^ {cot [(-1)^/^ y^ 



coth 



V_l)l/4 /z 



45 



(^)" 



45-^^ 



E 

n=l 



{0 



n^n^+i^ 



1 2 



(C2) 



(C3) 
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where Eqs. (|C2|) and (jC3|) arise from the contributions 
of |np and Rell^ , respectively. In the Hmit V jVc -^ oo , 
the sum of Eqs. (|C2|I and l)C3p saturates at 8/15 . Then 
there is the part of {SgSg)^ containing Re 11^ , which also 



does not contribute to the LK asymptotic behavior but 
gives a contribution which saturates towards a constant 
as the voltage is increased: 



Vl/Ve V2/Va 



{S9^9)2,Rc = - 16 / rfzi / dz2^TrRen^^ 

c 



- 32$: 



m^TT^ + m^n'^Tr'^ + n^-K^ + f -^ 



y— oo 



m,n>o ra^ ri^ 7i-4 



0.1905 




where again \n) — A/2sin(n7ry) . Finally, the |np contri- 
bution to {SgSg)^ will also be modified at finite V/Vc as 
compared to its asymptotics given by Eq. ((T7|l . The cor- 
responding expressionSi is too cumbersome, and we do 
not reproduce it here. Combining all the contributions 
we get the result shown graphically in Fig. ^ 



APPENDIX D: SUBLEADING TERMS IN THE 
STRONG-DEPHASING REGIME 

In this Appendix we show that the terms {Sgdg)^ 
and {Sg 6g) ^ remain small compared to the leading term, 
{Sg5g)2 , at high voltages, V/Vc S> g , i.e. in the regime 
of strong dephasing. Let us first analyze the behavior of 
the contribution H12|l . 



{SgSg),=32Vc f dz^ 



d „ 



^■„^Y2. 



V-,=V2 = V 



(Dl) 



which to leading order in Vc/V is 
d ~ 



{6g5g)^^32 



^v.-m 



Vl=V2 



dt E{t) 



(D2) 



The last term in the brackets is just twice {Sg Sg)^ , which 
we will discuss separately below. The first term in the 
brackets gives up to a numerical prefactor 



(S 9^9)1 



¥) 



5/4 



(D3) 



Clearly, Eq. IJD3|) is not only smaller by a factor 1/g 
compared to the leading term, Eq. (|48|l . at the lower 
bound of the strong dephasing regime V/Vc ~ 9 , but also 
decays faster with increasing voltage. 

It remains to show that the contribution of {Sgdg)^ 
also is small compared to Eq. H48f) . Setting Vi ~V2 = V 
in Eq. ((231, we get 



1 00 ^l(jl)=yi ?2(t2)=J/1 

Eo = 2 f dyidy2 f dTidT2 f V^i{ti) f P^fe) x 

ii{0)=y2 «2(0)=y2 

-Ti .2 T2 .2 min(Ti,r2) 

exp{- dti^ -~ dt2^ ~—F{y2)\Ti-T2\ -— rfi |6 - 6| 2/2(1 - y2) 

00' 



(D4) 



r 



Since the bare field fluctuations in the time interval be- 
tween Ti and T2 strongly suppress interference, relevant 
contributions are given by paths with ti and T2 close to 
each other. This allows us to approximate the path in- 



tegral by two trajectories propagating for the same time, 
keeping only the plain exponential suppression factor to 
account for the time difference: 
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?(ti)=0 



So = 2V2 fdyidy2 fdTJdr2 exp { -|^ ^(y2) In - raj 1 / V^{h) exp i - f dt 







C(o)=o 

r 






%/2iei 1/2(1 



V2) 



iM 



where again we transformed to the sum and the difference integral over the sum. We can also perform the integral 
of the coordinates ^1 , ^2 and carried out the trivial path over T2, 

I 



^0 



V2^ /• dy 



V J Fiy) 




dr 1 - e ^<^s 



1 00 



dr 



w-,^ya-y) 



1/3 



C(r)=0 

J Vat)expl-Jdt ^ + -^V2|e|y(l-y) 
e(o)=o 

T7 — J. I u ; / II 

1(0, r) 



1 — exp < 



V 



¥-any)r 



V^g 



w-g^yi^-y) 



2/3 



(D6) 



where {dr — d^ + |77|)X(77,t) — 6(1]) S{t) . This integral so that {Sg6g)Q remains smaller by a factor ~ 1/g than 
is dominated by y close to or 1, and should be cut off {6g Sg)r,_ in the strong dephasing regime. Since both 
at y ~ yc and {I — y) ^ y^ , where j/c is given by Eq. (|47|l . 
This yields 



{5g5g)^ 



9 



V 



Eq. (|D3p and Eq. (|D7p are small compared to Eq. (|1S|| at 
high voltages, we thus have shown that the asymptotic 
suppression of the conductance fluctuations is given by 

(D7) Eq. EHI)- 



Also at Petersburg Nuclear Physics Institute, 188350 

St. Petersburg, Russia. 
^ B. L. Altshuler, Pis'ma Zh. Eksp. Teor. Fiz. 41, 530 (1985) 

[JETP Lett. 41, 648 (1985)]. 
^ P. A. Lee and A. D. Stone, Phys. Rev. Lett. 55, 1622 

(1985). 
^ P. A. Lee, A. D. Stone and H. Fukuyama, Phys. Rev. B 

35, 1039 (1987). 
* A. L Larkin and D. E. Khmel'nitskii, Zh. Eksp. Teor. Fiz. 

91, 1815 (1986) [Sov. Phys. JETP 64, 1075 (1986)]; Phys. 

Scr. T14, 4 (1986). 
^ R. A. Webb, S. Washburn, and C. P. Umbach, Phys. Rev. 

B 37, 8455 (1988). 
^ S. B. Kaplan, Phys. Rev. B 38, 7558 (1988). 
'' D. C. Ralph, K. S. Ralls, and R. A. Buhrman, Phys. Rev. 

Lett. 70, 986 (1993). 
® R. Schafer, K. Hecker, and H. Megger, Phys. Rev. B 53, 

15964 (1996). 
® R. Haussler, H. B. Weber, and H. v. Lohneysen, J. Low 

Temp. Phys. 118, 467 (2000). 
^'^ C. Terrier, D. Babic, C. Strunk, T. Nussbaumer and 

C. Schonenberger, Europhys. Lett. 59(3), 437 (2002) and 

private communication. Experimental data discussed but 

not explicitly presented in this paper can be found in 



Ref.lU 

C. Terrier, Ph.D. thesis, Universitat Basel (2002) 
B. L. Altshuler and A. G. Aronov, in Electron- Electron In- 
teraction in Disordered Conductors, edited by A. L. Efros 
and M. PoUak (Elsevier, Amsterdam, 1985). 
B. L. Altshuler, A. G. Aronov, and D. E. Khmelnitsky, 
J. Phys. C 15, 7367 (1982). 

L. V. Keldysh, Zh. Eksp. Teor. Fiz. 47, 1515 (1964) 
[Sov. Phys. JETP 20, 1018 (1965)]. 

J. Rammer and H. Smith, Rev. Mod. Phys. 58, 323 (1986). 
Note that the diagrams c-diff and c-coop, as well as some 
possibilities of inserting additional impurity lines in the 
Hikami boxes, are missing in Fig. 1 of Ref. U However, 
the expression given in Ref. appears to give correctly the 
sum of all diagrams. 

Sh. Kogan, Electronic Noise and Fluctuations in Solids 
(Cambridge University Press, Cambridge, 1996). 
K. E. Nagaev, Phys. Rev. B 57, 4628 (1998). 
For a review, see Ya. M. Blanter and M. Biittiker, Phys. 
Reports 336, 1 (2000). 

L L. Aleiner and Ya. M. Blanter, Phys. Rev. B 65, 115317 
(2002). 
T. Ludwig, diploma thesis, Universitat Karlsruhe (2002). 



